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ABSTRACT 

We present contemporary infrared and optical spectroscopic observations of the 
type Iln SN 1998S covering the period between 3 and 127 days after discovery. 
During the first week the spectra are characterised by prominent broad H, He and 
C III/N III emission lines with narrow peaks superimposed on a very blue continuum 
(T~24000 K). In the following two weeks the C III/N III emission vanished, together 
with the broad emission components of the H and He lines. Broad, blueshifted ab- 
k>( i sorption components appeared in the spectra. The temperature of the continuum also 

. dropped to ~ 14000 K. By the end of the first month the spectrum comprised broad 

blueshifted absorptions in H, He, Si II, Fe II and Sc II. By day 44, broad emission com- 
ponents in H and He reappeared in the spectra. These persisted to as late as ~100-130 
days, becoming increasingly asymmetric. We agree with Leonard et al. (2000) that the 
broad emission lines indicate interaction between the ejecta and circumstellar mate- 
rial (CSM) emitted by the progenitor. We also agree that the progenitor of SN 1998S 
appears to have gone through at least two phases of mass loss, giving rise to two CSM 
zones. Examination of the spectra indicates that the inner zone extended to < 90 AU, 
while the outer CSM extended from 185 AU to over 1800 AU. 

We also present high resolution spectra obtained at 17 and 36 days. These spec- 
tra exhibit narrow P Cygni H I and He I lines superimposed on shallower, broader 
absorption components. Narrow lines of [N II], [O III], [Ne III] and [Fe III] are also 
seen. We attribute the narrow lines to recombination and heating following ionisation 
of the outer CSM shell by the UV/X-ray flash at shock breakout. Using these lines we 
show that the outer CSM had a velocity of 40-50 km/s. Assuming a constant veloc- 
ity, we can infer that the outer CSM wind commenced more than 170 years ago, and 
ceased about 20 years ago, while the inner CSM wind may have commenced less than 
9 years ago. During the era of the outer CSM wind the outflow from the progenitor 
was high - at least ~ 2 x 10~ 5 M Q yr . This corresponds to a mass loss of at least 
—0.003 Mq, suggesting a massive progenitor. The shallower, broader absorption is 
of width —350 km/s and may have arisen from a component of the outer CSM shell 
produced when the progenitor was going through a later blue supergiant phase. Al- 
ternatively, it may have been produced by the acceleration of the outer CSM by the 
© 0000 RAS radiation pressure of the UV precursor. 

We also describe and model first overtone emission in carbon monoxide observed 
in SN 1998S. We deduce a CO mass of - 10~ 3 M Q moving at -2200 km/s, and infer 
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1 INTRODUCTION 

Core-collapse supernovae (SNe) are believed to arise from 
massive progenitors (M>8-1OM0). The lower-mass (~ 8-10 
Mq) and higher- mass (> 20 Mq) progenitors experience 
heavy mass loss during the final stages of their evolution, 
several solar masses being ejected via a range of mass-loss 
rates. Consequently, in the vicinity of some core-collapse 
supernovae, dense circumstellar material (CSM) would be 
distributed according to the mass-loss history of the pro- 
genitor. The subsequent interaction of the freely-expanding 
supernova ejecta with the slowly-moving CSM generates a 
fast shock wave in the CSM and a reverse shock wave in 
the ejecta. The shocked regions emit high-energy radiation. 
The intensity of this emission depends on the density of 
the CSM and the ejecta, and the shock acceleration of the 
ejecta during the explosion (Chevalier & Fransson 1994). If 
the density of the CSM is small then the effects of the in- 
teraction only become significant several years after the ex- 
plosion when the supernova has faded. However, if the CSM 
near the supernova is relatively dense, strong CSM-ejecta 
interaction can begin shortly after the explosion. This is the 
case for type Iln (n=narrow line) supernovae. These events 
exhibit narrow emission lines in their spectra superimposed 
on broader emission profiles. They also exhibit a strong blue 
continuum (Schlegel 1990). However, the broad P Cygni ab- 
sorption components typical of normal type II SNe are weak 
or absent in type Iln SNe. The presence of variable, narrow 
line emission is a direct manifestation of the excitation of 
the dense CSM by the SN radiation. In addition, the pres- 
ence of broad Ha emission without a broad P Cygni ab- 
sorption component is a clear indication that the observed 
broad line emission is powered by the ejecta-wind interaction 
(Chugai 1990). Consequently, type Iln supernovae can pro- 
vide unique information about the progenitor and its later 
evolution, through the observed properties of their CSM. In 
addition, the study of the interaction of these supernovae 
with the CSM can provide vital clues about galaxy evolu- 
tion and the nature of active galactic nuclei (c.f. Terlevich 
et al. 1992). 

SN 1998S is the brightest type Iln event ever observed. 
It was discovered on 1998 March 2.68 UT in the highly- 
inclined Sc galaxy NGC 3877 by Z. Wan (Li & Wan 1998) 
at a broadband (unfiltered) optical magnitude of +15.2. The 
supernova is located at 16" west and 46" south of the nu- 
cleus. By March 18.4 it had brightened to V = +12.2 (Fassia 
et al. 2000). In a prediscovery frame obtained on 1998 Febru- 
ary 23.7 there was no evidence of the supernova, down to a 
limiting apparent magnitude of ~+18 (Leonard et al. 2000, 
IAUC 6835). We can thus assume that SN 1998S was discov- 
ered within a few days of the shock breakout. In the present 
paper we have adopted the discovery date 1998 March 2.68 
UT = JD 2450875.2 as epoch days, to, and express all 
other epochs relative to this fiducial date. 

Optical spectra obtained about day 3 (March 5-6) by 
Filippenko & Moran (1998), Huchra (Garnavich et al. 1998) 
and ourselves (see +3.3 d spectrum in Figure jij) showed 
prominent H and He emission lines with narrow peaks and 
broad wings superimposed on a blue continuum. As men- 
tioned above these narrow lines indicate the presence of a 
dense CSM in the vicinity of the supernova. Using opti- 
cal spectropolarimetry obtained at 5 days, Leonard et al. 



(2000) [LOO] deduced that the CSM is asymmetrically dis- 
tributed. Bowen et al. (2000) have presented high resolution 
spectroscopy of interstellar and circumstellar lines towards 
SN 1998S. They suggest that the CSM comprises a dense 
shell expanding at 50 km/s with a more highly ionised shell 
moving at ~300 km/s. They also estimate that SN 1998S is 
at a real distance of ~10 kpc from the nucleus and deduce 
from the interstellar lines that it lies on the far side of the 
galaxy disk. Gerardy et al. (2000) [GOO] presented near-IR 
spectra of SN 1998S spanning 95-355 days post-maximum 
light. They identified emission from carbon monoxide. They 
also suggested that late-time multi-peak H and He line pro- 
files in their optical and IR spectra indicate emission from 
a disk-shaped or ring-shaped circumstellar component. In 
addition, their t>225 d near-infrared spectra exhibit a con- 
tinuum that rises towards the longer wavelengths. They pro- 
pose that the rising continuum is likely due to dust heated 
by the interaction of the ejecta with the CSM. 

In Fassia et al. (2000) we presented contemporaneous 
optical and infrared photometric observations of SN 1998S 
covering the period between 11 and 146 days after discovery. 
Using the interstellar Na I D lines we derived an extinction 
Av = 0.681q 25 ma g- We also examined the evolution of 
the total luminosity and found that during the first month 
the luminosity decreased very rapidly viz. ~0.08 mag/day. 
Subsequently (30-70 days) the decline rate decreased, re- 
sembling that of type IIL supernovae viz. 0.05 mag/day. By 
day ~100 the decline rate slowed to 0.01 mag/day matching 
the radioactive decay of 56 Co. From the bolometric lumi- 
nosity after ~100 we estimate that 0.15+0.05 M of 56 Ni 
were produced in the explosion. Furthermore, we discovered 
that as early as day 130 the supernova exhibited an aston- 
ishingly high infrared (IR) excess, K — L' = +2.5. We argue 
that this excess is due to dust grains in the vicinity of the 
supernova. However, the physical extent of this early IR lu- 
minosity source was so large that the emission must have 
come from pre-existing dust in the CSM, possibly powered 
by X-rays from the ejecta-CSM interaction. 

In this paper we present optical and infrared spec- 
troscopy of SN 1998S covering the period 3-127 days after 
discovery. The observations and the reduction procedure are 
presented in section 2. In section 3 we discuss possible line 
identifications. In section 4 we present an overview of the 
spectral behaviour of SN1998S and derive constraints about 
the nature and the characteristics of the CSM. We also anal- 
yse and model the first overtone of the CO emission. The 
work is summarised in section 5. 



2 OBSERVATIONS 

2.1 Optical spectroscopy 

2.1.1 Low resolution spectroscopy 

Optical spectroscopy of SN 1998S was acquired using the 
ISIS spectrograph and the UES echelle spectrograph at 
the William Herschel Telescope (WHT), La Palma, the 
IDS spectrograph at the Isaac Newton Telescope (INT), 
La Palma and the fiber-fed HYDRA spectrograph at 
the Wisconsin-Indiana- Yale-NOAO Telescope (WIYN), Kitt 
Peak. The log of optical observations is shown in Table 

Apart from the 17.4, 36.3d UES spectra, the spectro- 
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Table 1. Log of optical spectroscopy of SN1998S 
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a 2450000+ 

b Used to correct for atmospheric and instrumental transmission and to flux calibrate the supernova spectra. 
c This is a fibre-fed spectrograph. The diameter of the fiber used was 3.1 arcsec. 
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Figure 1. Optical spectra of SN 1998S taken at the Isaac Newton Telescope (La Palma) and the Wisconsin- Indiana- Yale-NOAO Telescope 
(Kitt Peak) (see Table [l] for details). The data have not been corrected for redshift or reddening. The epochs are with respect to the 
discovery date (1998 March 2.68 UT= days). For clarity, the spectra have been displaced vertically by the amounts indicated (in units 
of 10 -14 ergs s — 1 cm -2 A -1 ). The 97. 2d spectrum has also been multiplied by 2. 
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Figure 2. High-resolution (echelle) optical spectra of SN 1998S taken at the William Herschel Telescope (WHT) on days 17.4 and 36.3. 
Also shown as a dotted line on the top-right plot is the high-resolution spectrum of the H a line obtained at the Wisconsin-Indiana- Yale- 
NOAO (WIYN) Telescope on day 38.5(see Table Q for details). The spectra have not been corrected for redshift. The Na I A5889, 5895 
night sky lines are marked with an asterisk. The vertical dashed line corresponds to the adopted redshift for the SN centre of mass, 
847±3 km/s (cf. Section 3.2) 
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scopic data were reduced by means of standard routines 
in the data-reduction package FIGARO (Shortridge 1995). 
The data were first cleaned of cosmic rays. Subsequently, 
the frames were debiased, fiat-fielded and sky-subtracted. 
The spectra were extracted from the resulting frames using 
simple extraction. Wavelength calibration was performed by 
means of Cu-Ar and Cu-Ne lamps. The spectrophotometric 
flux standards used are listed in Table jj]. 

Owing to adverse weather conditions and the use of 
fairly narrow slits (slit width < 1.6") in most cases, the 
fluxing described above was not accurate. Therefore to im- 
prove the fluxing, we made use of our BVR photometry of 
SN 1998S (Fassia et al. 2000) to derive correction scaling 
factors. Transmission functions for the B, V and R bands 
were constructed using the JKT filter functions, the CCD 
response and the standard La Palma atmospheric transmis- 
sion. The spectra were then multiplied by these functions 
and the resulting total flux within each band compared with 
the observed magnitudes. For dates where simultaneous pho- 
tometric observations did not exist we interpolated within 
the photometric data (Fassia et al. 2000). However, we do 
not have any photometry prior to 11 days. Therefore, to 
scale the 3.3 d spectrum, we used BVR photometry (day 5) 
from Garnavich et al. (1999) together with U-band and unfil- 
tered CCD photometry from the IAU circulars (IAUC 6829, 
IAUC 6831, IAUC 6835) covering days 0-7.5 to interpolate 
to day 3.3. For any particular epoch the scaling factors de- 
rived for each band agreed to within ±10% demonstrating 
good internal consistency for the procedure. The spectra 
were multiplied by the geometric mean of the scaling fac- 
tors for each epoch. The amounts by which the spectra had 
to be scaled to match the photometry ranged from xl.23 to 
X2.48. Including uncertainties in the photometry, we esti- 
mate the final absolute fluxing of the spectra to be accurate 
to <15%. The optical spectra are shown in Figure [l]. 

2.1.2 High resolution spectroscopy 

High resolution optical spectra were obtained on days 17.4 
and 36.3 with the WHT UES echelle spectrograph, and on 
day 38.5 with the WIYN Hydra spectrograph using a high 
dispersion grating. Some of the UES spectra have already 
been presented and discussed in Bowen et al. (2000). The 
WHT/UES high- resolution spectra were reduced using the 
IRAF package ECHELLE. The data were bias-subtracted 
and flat-field corrected. The different echelle orders were 
then traced and extracted from the frames. Wavelength 
calibration was by means of a Thorium- Argon lamp. We 
checked the precision of this calibration using night-sky 
emission lines present in our spectra together with the high- 
resolution night-sky line atlas of Osterbrock et al. (1996). 
This demonstrated that our wavelength calibration is accu- 
rate to ±2 km/s. 

Cosmic rays were present in the extracted spectra. 
These were identified and removed by comparing repeat ob- 
servations of the spectra obtained at different times during 
the night. The counts in the contaminated pixels were then 
replaced by those of the corresponding pixels in the other 
spectrum. The mean fluxes of the spectra were then nor- 
malised to unity by fitting a low-order spline function to the 
instrumental continuum. Finally the spectra were fluxed as 
follows. First we determined the shape and flux of the over- 



all SN continuum in the optical region at days 17.4 and 36.3. 
This was done by fitting a low-order spline function to the 
continuum in contemporary low-resolution spectra. This was 
straightforward for the earlier epoch since a low-resolution 
spectrum was available on day 17.7. However, for day 36.3 
no contemporary low-resolution spectrum was available. We 
therefore used the 27.2 d low-resolution spectrum scaled to 
match the BVR magnitudes on day 36.3d. Given the evolu- 
tion of the B — V and V — R colours, we can argue that the 
continuum shape did not change much during this period. 
We then scaled the normalised high-resolution spectra so 
that the continuum in the regions where lines were observed 
matched the continua determined from the low-resolution 
spectra. We estimate the fluxing of the high resolution lines 
to be accurate to better than ±20%. In Figure || we show 
the spectral regions where lines are observed. Also shown in 
Figure | is the high-resolution 38. 5d WIYN /HYDRA spec- 
trum. The HYDRA spectra were reduced in the same way as 
the low- resolution optical spectra (cf Section 2.1) and were 
fluxed following the same procedure as for the day 36.3 high- 
resolution spectra. 

2.2 Infrared Spectroscopy 

Infrared spectra of SN 1998S were obtained at the United 
Kingdom Infrared Telescope (UKIRT), Hawaii, using the 
Cooled Grating Spectrograph CGS4. The observing log is 
shown in Table iB. During the observations the telescope was 
nodded 7.5 arc seconds along the slit. The data were re- 
duced using CGS4DR (Daly & Beard 1992) and FIGARO. 
After bias-subtraction and flat-field correction the pair of 
images obtained at the two nod positions were subtracted 
to remove sky line emission. If residual sky line emission 
persisted, it was removed using POLYSKY. The positive 
and negative spectra were then extracted from the result- 
ing frames using the optimal extraction algorithm of Home 
(1986). Wavelength calibration was by means of argon and 
krypton arc lamps. Observations of standard stars were used 
to correct for the atmospheric and instrumental transmis- 
sion and to flux-calibrate the spectra. Hydrogen absorption 
lines in the standard star spectra were removed before these 
corrections/calibrations were performed. The standard stars 
used are listed in Table |i| 

As with some of the optical spectra, variable weather 
conditions together with the use of narrow slits for some 
of the spectra meant that fluxing using standard stars was 
approximate only. To improve the fluxing we could, in prin- 
ciple, follow a similar procedure to that applied to the opti- 
cal spectra. However, no IR photometry was available until 
13.8 d, 11 days after our earliest IR spectrum. Moreover, for 
most of the J band spectra the J band photometry trans- 
mission function (filter±atmosphere±detector response) ex- 
tended significantly beyond the red limit of the spectra. 

For the earliest IR spectrum, which was at 2.8 d in 
K, there is a particular problem in checking the fluxing. 
Not only is there a lack of accurate photometry around this 
time, but also optical observations indicate that the light 
curve was rising rapidly, making it very difficult to make a 
flux correction using extrapolated magnitude values. How- 
ever, the spectrum was taken in a wide (2.46") slit and so 
fluxing by means of the standard star is accurate to ±20%. 
Therefore, no attempt was made to improve this fluxing. 
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Table 2. Log of infrared spectroscopy of SN1998S obtained at UKIRT with CGS4 
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a 2450000+ 

b Used to correct for atmospheric and instrumental transmission and to flux calibrate the supernova spectra 
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Figure 3. Blackbody fits (t <44.6 days) and spline fits (t >44.6 days) to UBVRIJHK photometry (see section 2.2 for details). These 
fits were used to ascertain the absolute fluxing of the infrared spectra. 
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Figure 4. J band spectra of SN 1998S taken at the United Kingdom Infrared Telescope (UKIRT). The spectra have not been corrected 
for redshift. The epochs are with respect to the discovery date (1998 March 2.68 UT= days). For clarity, the spectra have been displaced 
vertically by the amounts indicated (in units of 10 -15 ergs s _1 cm~ 2 A -1 ). 
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Figure 5. H and K band spectra of SN 1998S taken at the United Kingdom Infrared Telescope (UKIRT). The spectra have not been 
corrected for redshift. The epochs are with respect to the discovery date (1998 March 2.68 UT= days). For clarity, the spectra have 
been displaced vertically by the amounts indicated (in units of 10 -15 ergs s _1 cm~ 2 A -1 ). 
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For the other IR spectra we corrected the fluxing as 
follows. Using the optical and IR magnitudes given in Fas- 
sia et al. (2000), together with earlier optical magnitudes 
reported by Garnavich et al. (1999) we interpolated to pro- 
vide magnitudes for each epoch upon which IR spectra were 
acquired. Thus, a set of UBVRI magnitudes were estab- 
lished for epochs 4.8 d to 8.9 d, and UBVRIJHK magni- 
tudes for epoch 11.9 d onward. The magnitudes were then 
converted to fluxes using the calibrations of Wilson et al. 
(1972) and Bessel (1979). For epochs 4.8 d to 8.9 d we fit- 
ted reddened blackbody functions to the UBVRI data. The 
blackbody functions were reddened using the empirical for- 
mula of Cardelli et al. (1989) with A v = 0.68i^2B (Fassia et 
al. 2000). We then scaled the fluxes of the infrared spectra 
so that the continua matched the blackbody functions (cf. 
Figure |^) . This procedure ignores the contribution of the op- 
tical emission lines to the blackbody fits, probably leading 
to an overestimate of the IR flux. However, from the 3.3 d 
optical spectrum, we estimate this systematic error was no 
more than 7%. For epochs 11.9 d to 17.7 d, while contempo- 
raneous IR photometry was available, correction of the J- 
band spectra was still difficult due to the incomplete overlap 
problem explained above. Therefore, for the 11.9 and 17.7 d 
J-band spectra we followed a similar procedure to that car- 
ried out for the earlier epochs, except that now we included 
JHK photometry in the blackbody fits. Given the slow vari- 
ation of the supernova flux around this time {cf. Fassia et al. 
2000, Garnavich et al.), we simply scaled the 11.9 d J spec- 
trum so that its continuum matched the reddened blackbody 
14.3 d optical spectrum. The 12.8 d 18834-25240 A spectrum 
was also flux-corrected by matching it to the same reddened 
blackbody. This spectrum was also checked using contem- 
poraneous H and K light curves (see below). Similarly, the 
17.7 d J-band spectrum was matched to the 17.7 d optical 
spectrum. 

By epoch 44.6 d and later, the IR spectra were increas- 
ingly affected by emission/absorption features (cf. Figures [jj 
^| [|). Consequently, flux correction using the blackbody fit- 
ting procedure was increasingly inappropriate. Instead, for 
these later epochs we applied spline fits (Press et al. 1992) 
to the UBVRIJHK photometry. We then scaled the in- 
frared spectra so that the integrated fluxes of the observed 
IR spectra and spline spectra were equal in the wavelength 
regions where they coincided (cf. Figures ^|). Scaling factors 
for the infrared spectra ranged from xO.92 to xl.65. 

To check the above scaling procedures we used the op- 
tical flux correction method (Section 2.1) for the few IR 
cases where this was possible. Thus, using contemporary IR 
photometry (Fassia et al. 2000), we estimated absolute flux 
levels for the 12.8 d and 108. 6d K spectra and the 126.6 d 
H spectrum. We constructed transmission functions for the 
H and K bands using the IRTF/UKIRT filter functions 
and the standard Mauna Kea atmospheric transmission. We 
then multiplied the corresponding spectra by these functions 
and compared the resulting total fluxes with the photome- 
try. The scaling factors derived by this procedure agreed to 
within 12% with those derived by the blackbody /spline fits 
mentioned above. We estimate, therefore, that the absolute 
fluxing of the corrected infrared spectra is accurate to <15%. 
The final fluxed infrared spectra are shown in Figures KM and 
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On days 44.6 and 55.6, high resolution (~100 km/s) 
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Figure 6. High-resolution spectra of SN 1998S taken at the 
United Kingdom Infrared Telescope (Hawaii). The epochs are 
with respect to the discovery date (1998 March 2.68 UT= days). 
For clarity, the spectra have been displaced vertically by the 
amounts indicated (in units of 10 — 15 ergs s _1 cm -2 A -1 ). 



spectra of the He I 10830 A line were acquired with 
UKIRT/CGS4. They were wavelength calibrated using sky 
lines from the raw images. Fluxing was carried out by scal- 
ing their continua to match the contemporary low-resolution 
spectra. They are shown in Figure ^. 



3 DESCRIPTION OF THE SPECTRA 
3.1 Low-Resolution Spectra 

In Figures [j], ^| we redisplay the optical and IR spectra, re- 
spectively, but with the best blackbody/spline fit continua 
(see Section 2.2) subtracted in order to minimise distortion 
of the line profiles. In addition the relative flux is plotted 
logarithmically in order to bring out the detail in the line 
profiles. The spectra have not been redshift corrected. The 
resulting line identifications are indicated in Tables 3, 4 and 
in Figures (?], [| In general, the lines have complex profiles 
comprising mixtures of broad and narrow emission and ab- 
sorption lines. To examine these profiles in detail we car- 
ried out multiple component fits using Gaussian profiles. 
This provided peak wavelengths, intensities and widths of 
the separate components of each profile. These are listed in 
Tables 3, 4 for the optical and infrared spectra respectively. 
We note that the presentations of LOO and GOO include first 
season optical and IR spectra which have similar epochs to 
some of the spectra discussed here. 



3.1.1 3-7 days 

During the first week, the spectra are characterised by 
strong, complex emission lines of the Balmer, Paschen 
and Brackett series, He I 10830 A, He II 4686 A and 
C III 4648/N III 4640 A, lying on a smooth blue continuum. 
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There is little sign of any absorption components. In Fig- 
ure | we show blackbody spectra reddened by E(B-V)=0.22 
(Ry=3.1) superimposed on the optical and IR continua. The 
blackbody temperature on day 3.3 is 23,600±500 K. A some- 
what higher blackbody temperature of 28,000 K is derived 
by LOO for their day 4 optical spectrum, using a slightly 
higher E(B-V) of 0.23. They point out that this is quite 
high for an SN II at this epoch. In agreement with LOO, we 
also find that the day 3-4 blackbody fit fails to account for 
all the continuum blueward of ~5000 A. LOO argue that the 
most likely cause is a departure from the blackbody spec- 
trum due to the increase of the continuum absorptive opac- 
ity (bound-free, free-free) with wavelength. We confirm the 
LOO finding that a typical line profile at this time comprises 
a narrow unresolved component superimposed on a broad 
base (cf. Figures ^). This is illustrated in Figures ^| and 
|ic| where we show multiple gaussian component fits to Ha 
and He I 10830 A respectively. The unresolved component 
has a FWHM < 500 km/s, superimposed on a broad compo- 
nent of FWHM ~4000 km/s. Other broad emission features 
include one at ~7100 A which we attribute to a blend of 
He I 7065 A and C II 7100 A. Other narrow lines include 
C III 5696 A, He I 5876 A, 6678 A, 7065 A and He II 5411 A, 
8236 A (possibly) and 11626 A. Note that the narrow ab- 
sorption feature at 20,000 A is residual absorption due to 
H2O, CO2 and CH4 vapour in the earth's atmosphere. 



3.1.2 8-18 days 

Between days 3.3 and 17.7, blackbody fits indicate a cool- 
ing continuum. By day 17.7 the blackbody temperature was 
14500 K, consistent with the 13,400±1200 K derived by 
Fassia et al. (2000) from contemporary photometry. Dur- 
ing this time, the lines also underwent a substantial change. 
Between days 3.3 and 14.3 the broad emission component 
in Balmer, Paschen, Brackett and He I 10830 A lines largely 
vanished (see also Qiu et al. 1998), leaving behind nar- 
rower, but still quite complex, emission features. The intense 
C III 4648/N III 4640 A feature disappeared altogether. 
The disappearance of the broad emission in the Balmer and 
He I 10830 A lines was accompanied by the appearance of 
a broad, blueshifted absorption extending to velocities of 
~10000km/s. However, broad absorptions did not appear 
in the Paschen or Brackett lines. In Ha and He I 10830 A 
weak, unresolved (FWHM < 500 km/s) P Cygni profiles de- 
veloped. These narrow lines are examined in greater detail 
in section 3.2. 



3.1.3 27-56 days 

The day 27 optical spectrum is very similar in appearance 
to the day 25 spectrum of LOO. In making line identifica- 
tions, we were guided by those given for the photospheric 
spectra of the type II SN 1987A and SN 1995V (Williams 
1987, Fassia el al. 1998). By day 27 the spectrum com- 
prises broad, blueshifted absorptions in H (Balmer series), 
He I+Na I blend, Si II, Fe II and Sc II, together with narrow 
P Cygni lines of H (Balmer series), He I and narrow emis- 
sion lines of [O III] . (Narrow emission lines of [N II] , [Ne III] 
and [Fe III] were also seen around this time in our high res- 
olution spectra - see section 3.2). The FWHM of the broad 



absorptions is ~4000 km/s, with the minima blueshifted by 
4000-5500km/s with respect to the local standard of rest. 
The blue edge of the Ha absorption trough is blueshifted by 
~8000 km/s. In the IR spectra, by days 44 and 55, broad 
emission features (FWHM ~5000km/s) had developed in 
the Paschen, Brackett and He I 10830 A lines. The FWHM 
of the Paa, Pa/3 and Br7 lines on day 44.6 is ~5000 km/s. A 
broad, blueshifted absorption in He I 10830 A is also appar- 
ent. The deep narrow absorption seen in He I 10830 A on 
day 17 is still present on days 44 and 55. This is discussed 
in section 3.2. Other broad IR features that had appeared 
by this period we identify with blends of O I 11287/94- 
S I 11306 A, and Mg I 11828-Si II 12047 A (cf. Meikle 
el al. 1989). LOO note that the relative strengths of some 
of the broad features in the optical spectrum resemble those 
found in the spectra of SNe Ic, suggesting hydrogen-deficient 
ejecta. They also suggest that the low contrast of the broad 
lines relative to the continuum may be due to external illu- 
mination of the line-forming region by light from the CSM- 
ejecta interaction (see also Lentz et al. 2000 and Branch et 
al. 2000). This is discussed in more detail in Section 4.1.2. 

3.1.4 72-127 days 

Between days 27 and 72 in the optical region the most strik- 
ing change is the replacement of the fairly weak, broad 
P Cygni Ha feature with a strong, broad (FWHM ~7000 
km/s) emission feature (cf. Figure [j]). The profile is some- 
what asymmetric. Its peak is blueshifted by ~-500km/s with 
respect to the SN rest frame, and its blue wing is steeper. 
The extreme limit of the blue wing is blueshifted by ~-8500 
km/s. The observation of comparably wide Paschen emis- 
sion lines on day 44.6 suggests that broad Ha emission had 
probably already appeared by then. In other Balmer lines, 
the change from days 27 to 72 was less dramatic. The ab- 
solute depth of the broad absorption in H/3 weakened and 
the H7 absorption disappeared almost entirely by day 97. 
A dramatic change between days 27 and 72 is the appear- 
ance of strong, broad emission in the Ca II IR triplet. The 
presence of O I 7773 A and O I 9264 A emission lines on 
the 72 d spectrum suggests that significant emission from 
the O I 8446 A line may be blended with the Ca II mul- 
tiplet. The broad absorption lines in the He I+Na I blend, 
Fe II and Sc II became more pronounced, but vanished in 
Si II 6355 A. The narrow H I and He I lines had weakened 
by day 72, with only Ha being clearly visible. Weak narrow 
[O III]4959, 5007 A lines are visible in the day 97 spectrum. 

Our day 97 optical spectrum is similar to the day 108 
spectrum of LOO, and our day 109/113/127 composite IR 
spectrum is similar to the day 110 spectrum described in 
GOO. By this period, the spectra were dominated by broad 
emission lines of Ha, Paa, Pa/3, Pa7+He I 10830 A blend 
and the Ca II triplet. The Ha, Paa and Pa/3 lines had 
broadened to a FWHM of 7000-8500 km/s. Moreover the 
Ha, and possibly Paa, profiles had become more asymmet- 
ric, with the development of a pronounced steep blue edge. 
The asymmetry, and the pronounced red wings of these lines 
could indicate that Thomson scattering of line radiation oc- 
curred in the ejecta during this epoch. We note however 
that Br7 had weakened significantly. Also of importance is 
the appearance by day 109 of first overtone emission from 
carbon monoxide. This is discussed further in section 4.2. 
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Figure 7. Optical spectra of SN 1998S with blackbody/spline fit continua shown in Figure |i] subtracted. Suggested line identifications 
are indicated. (The dashed lines correspond to the rest wavelength of the lines in NGC 3877). Residual telluric absorption lines are marked 
with the dotted lines and an asterisk. 



LOO suggest that the appearance of broad, asymmetric Ha 
emission in their day 108 spectrum is due to interaction of 
the ejecta with that part of the CSM responsible for the nar- 
row lines seen at earlier epochs. We support this proposal, 
and suggest that the presence of strong, broad Paschen and 
Brackett lines on day 44 indicates that the interaction may 
have already begun as early as this. This will be discussed 
further in Section 4. 



3.2 High-Resolution Spectra 

The WHT/UES observations of SN 1998S on days 17.4 and 
36.3 revealed narrow absorption and emission lines. Many 
of the absorption lines have widths of around 2-15 km/s 
and can be attributed to the interstellar medium (ISM) of 
the Milky Way and NGC 3877. Fassia et al. (2000) used the 
Na I D absorption lines observed on day 36.3 to estimate the 
extinction to SN 1998S. In a more comprehensive analysis 
of the WHT/UES ISM data, Bowen et al. (2000) examined 
a range of lines in the UV and optical regions. They also 
describe a blueshifted absorption feature in Mg II 2796 A 
and in the Balmer lines, of width ~350 km/s. Variation in 
the depth of these features between different epochs points 
to an origin in the CSM. A variable, narrow component in 
Na I D is also attributed to the CSM. In addition, Bowen et 
al. report the presence of narrower (~30 km/s) Balmer and 
optical He I P Cygni lines as well as similarly narrow P Cygni 
lines in a number of Fe II UV multiplets. These too are 
attributed to the CSM of SN 1998S. They also indicate (in 



their Table 4) the presence on day 36.3 of narrow forbidden 
lines of [N II], [Ne III], [O III] and [Fe III], and suggest an 
origin in either the CSM or some other line-of-sight emission 
line region. 

In the present work, we consider only those narrow lines 
whose origin is in the CSM of SN 1998S. In Figure | we dis- 
play and identify all the CSM lines present in the day 17.4 
and 36.3 WHT/UES spectra (resolution ~7 km/s). Also 
shown is the day 38 Ha profile (resolution ~12 km/s) ob- 
tained at the WIYN. In Figure]^ we show the high resolution 
(~100 km/s) P Cygni spectra of the He I 10830 A line ac- 
quired with UKIRT/CGS4 on days 44.6 and 55.6. In the 
optical spectra we confirm the presence of most of the al- 
lowed and forbidden lines described in Bowen et al. (2000). 
In particular, we confirm that the H lines are composed of 
two components viz. a narrow P Cygni-like profile superim- 
posed on a broader absorption feature which extends as far 
as ~-350 km/s to the blue. We also note the presence of 
a similar broad absorption in He I 5876 A and possibly, in 
He I 4471, 5016 A. Identifications of the narrow lines are 
listed in Table |B| along with the redshifts of the emission 
peaks, the blueshifts of the P Cygni troughs relative to the 
estimated SN centre of mass (see below), the widths of the 
emission and absorption components and the line intensi- 
ties. The narrow H profiles span a range of ~130 km/s, the 
He I profiles span ~100 km/s, while those of the forbidden 
lines span ~80km/s. 

Two of the narrow emission lines which Bowen et al. 
(2000) list in their Table 4, they identify as He I 3870 A 
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Figure 8. Infrared spectra of SN 1998S with the blackbody /spline fit continua shown in Figure ^subtracted. Suggested line identifications 
are indicated. (The dashed lines correspond to the rest wavelength of the lines in NGC 3877) 



and O I 7773 A. We confirm the presence of these features. 
However, it seems unlikely that they originate in the SN 
CSM. Firstly, both lines are unresolved, implying a linewidth 
of less than ~10 km/s. This is much narrower than the 
FWHM~60 km/s which is typical of the emission compo- 
nents of other narrow lines. The "He I 3870 A" feature can 
be seen in Figure | lying ~100 km/s redward of [Ne III] 
3869 A. In the case of this line, comparison of two succes- 
sive repeat frames shows that it appeared strongly only in 
the first, while in the second it was not discernible above 
the noise. In the case of the "O I 7773 A" line we note that 
it extends the full length of the slit (15") equivalent to over 
1 kpc at NGC 3877. It is also at the wrong redshift for it 
to be associated with SN 1998S. It therefore seems unlikely 
that these two lines originate in the SN CSM. We also note 
that in Table 4 of Bowen et al, "He I 4290 A" should proba- 
bly be read as "He I 4920 A" . Bowen et al. also list a narrow 
P Cygni line identified with He I 4713 A. While there is a 
weak feature at this wavelength in our reduced data, its S/N 
is very low and so we do not include it in Table ^| Narrow 
optical lines identified by us but not listed by Bowen et al. 
include He I 6678, 7065 A and [Fe III] 4658 A. 

The forbidden lines on day 36.3 all exhibit about 
the same absolute redshift and width (see below), with a 
weighted mean redshift of +846.9±2.9 km/s. We adopt this 
as probably being closest to the SN centre of mass red- 
shift, for the following reasons. We note below that between 
days 17.4 and 36.3 the red wings of the emission components 
of the H lines (and possibly also the He I lines) and forbidden 
lines moved to the red by ~+ 10-20 km/s. We suggest that 



this is a manifestation of light travel time effects across the 
CSM around the supernova. The radius of the outer CSM is 
at least 10 light days (see below) and so on the timescale of 
the observations as the SN-flash paraboloid moves further 
into the CSM on the far side of the SN, increasing amounts 
of more redshifted material will contribute to the emission 
lines. We therefore favour the later epoch (36.3 days) since 
by that time the emission lines should be originating from a 
greater proportion of the CSM i. e. the bias to the blue due 
to light travel time effects should be reduced. We also favour 
the forbidden lines since the allowed lines appear to be con- 
siderably more complex, comprising a mixture of broad and 
narrow P Cygni components plus a probable additional re- 
combination emission component. Our adopted redshift for 
the supernova CSM, +847 km/s, is consistent with the ap- 
proximate value of +860 km/s indicated in Bowen et al. 
(2000). We note that measurements of the systemic veloc- 
ity of NGC 3877 vary from +838+11 km/s using optical 
observations of the nuclear region (Fouque et al. 1992), to 
902+6 km/s (De Vaucouleurs et al. 1991) and +910+6 km/s 
(Broeils & van Woerden 1994) using an H I 21 cm emission 
line map. 

We have examined variations in shifts, widths and inten- 
sities of the stronger narrow lines between different species 
and epochs. The lines considered were Hq,/3,7, He I 4471, 
5015, 5875 A, [O III] 4959, 5007 A, [Nil] 5755 A, [Ne III] 
3869 A and [Fe III] 4658, 5270 A. In Figure [H] we show 
in detail those lines which were detected on both days 17.4 
and 36.3 viz. H, He I and [O III]. The relative velocity shifts 
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FLUX' 


He 


BE 








4001.2(2.6) 


28.2(4.1) 


0.80(0.25) 




























NE 








3985.3(0.3) 


4.4(0.7) 


0.15(0.03) 




























NA 








3978.1(0.3) 


4.5(0.3) 


0.22(0.02) 


























Sr II 4077 


NA 




















4074.8(1.4) 


8.8(1.9) 


0.12(0.04) 














H5 


BE 


4108.7(5.1) 


84.1(7.3) 


4.02(1.10) 


4141.8(3.3) 


44.1(8.6) 


0.75(0.35) 




























BA 




















4072.0(7.9) 


93.3(9.4) 


1.75(0.35) 
















NE 


4108.8(2.0) 


24.1(4.5) 


1.16(0.37) 


4115.6(0.4) 


8.6(0.7) 


0.27(0.02) 




























NA 








4097.0(0.4) 


8.7(1.5) 


0.20(0.05) 








4107.2(2.1) 


13.3(3.9) 


0.25(0.11) 














Hy 


BE 


4353.8(4.9) 


69.9(13.7) 


2.68(0.57) 


4370.3(4.7) 


45.6(9.7) 


1.17(0.35) 


4365.9(2.8) 


32.7(11.2) 


1.54(0.35) 






















NE 


4349.5(2.0) 


12.0(4.9) 


0.6(0.2) 


4354.3(0.3) 


6.9(0.5) 


0.25(0.03) 




























BA 




















4271.4(5.2) 


103.0(14.5) 


5.96(1.9) 
















NA 








4343.3(2.0) 


19.2(2.9) 


0.29(0.07) 








4344.3(1.5) 


16.8(3.2) 


0.38(0.14) 














He 1 4471 


BA 




















4468.6(11.2) 


130.2(26.0) 


8.0(1.7) 














He II 4686 


BE' 


4677.8(2.6) 


101.7(5.6) 


15.2(1.06) 


































NE 


4696.3(1.2) 


14.0(1.2) 


1.25(0.16) 
































C III 4648/N III 4640 


NE 


4650.9(0.6) 


15.9(0.5) 


1.4(0.08) 
































HP 


BE 


4873.0(2.2) 


62.6(7.5) 


3.98(0.34) 


4888.5(1.7) 


26.4(1.8) 


0.65(0.09) 


4893.7(10.4) 


33.2(8.2) 


0.67(0.10) 






















NE 


4873.4(0.5) 


9.2(0.9) 


1.02(0.11) 


4875.6(0.2) 


9.4(0.4) 


0.44(0.03) 


4876.6(1.7) 


10.6(4.9) 


0.35(0.18) 


4882.8(4.7) 


17.8(4.9) 


0.30(0.16) 
















BA 




















4802.7(11.2) 


95.7(18.1) 


4.23(1.9) 


4765.8(6.1) 


122.2(15.6) 


2.05(0.25) 


4753.9(4.5) 


64.5(15.5) 


0.31(0.08) 




NA 








4867.1(0.7) 


9.2(1.1) 


0.17(0.04) 


4869.6(1.2) 


6.0(1.2) 


0.22(0.07) 


4870.4(1.7) 


9.9(0.9) 


0.39(0.08) 














l ; e II 5018 


BA 




















4965.6(11.2) 


96.3(14.6) 


3.06(0.7) 


4974.4(1.1) 


56.6(14.5) 


1.38(0.15) 


4982.4(2.6) 


45.8(10.2) 


0.21(0.09) 


Fell 5169 


BA 




















5106.1(6.0) 


70.1(10.7) 


1.60(0.32) 


5118.1(1.3) 


80.8(3.2) 


2.17(0.5) 


5123.9(4.6) 


86.9(11.5) 


0.56(0.18) 


Hell 5411 


NE 


5419.7(1.2) 


23.7(4.0) 


0.38(0.07) 
































Sc II 5527 


BA 




















5434.9(5.5) 


79.7(14.5) 


0.48(0.90) 


5477.2(3.6) 


71.7(17.4) 


0.85(0.12) 


5474.2(27.8) 


181.0(50.7) 


0.53(0.19) 


Sc II 5658 


BA 


























5624.6(4.4) 


87.4(28.1) 


1.74(0.40) 


5638.8(6.8) 


111.4(11.6) 


0.63(0.16) 


CIII 5696 


NE 


5710.1(1.2) 


12.5(1.3) 


0.21(0.03) 
































He 1 5876 


BE* 


























5913.4(7.6) 


92.7(20.3) 


1.73(0.60) 










NE 


5892.1(1.2) 


8.50(0.8) 


0.16(0.02) 


5885.4(0.6) 


10.2(1.2) 


0.19(0.08) 


5894.5(0.8) 


8.8(1.8) 


0.22(0.08) 






















BA' 




















5806.5(3.1) 


91.3(9.1) 


2.53(0.21) 


5823.1(10.5) 


115.2(30.2) 


3.84(1.30) 


5807.9(3.5) 


133.9(24.6) 


1.14(0.30) 




NA 








5894.5(0.4) 


9.9(1.7) 


0.25(0.06) 


5886.9(1.2) 


8.3(2.7) 


0.13(0.04) 


5886.1(1.3) 


11.8(1.9) 


0.28(0.06) 














r> 


BA 


























6081.9(5.2) 


93.8(12.5) 


0.64(0.16) 


6065.5(3.0) 


148.4(12.5) 


0.64(0.14) 


Sc II 6245 


BA 


























6201.7(2.8) 


51.9(12.4) 


0.36(0.11) 


6210.3(5.0) 


68.7(16.8) 


0.16(0.06) 


Si II 6355 


BA 




















6286.8(2.3) 


81.7(11.6) 


2.49(0.51) 














Ha 


BE 


6579.4(5.2) 


97.5(15.7) 


5.7(0.9) 


6592.9(1.7) 


27.1(8.2) 


0.54(0.17) 


6599.5(10.0) 


37.7(15.5) 


0.91(0.35) 








6571.1(5.3) 


149.0(9.1) 


17.42(0.72) 


6581.7(3.2)/6602.6(3.2/ 




11.27(0.90) 




NE 


6581.6(2.6) 


11.9(1.9) 


1.44(0.22) 


6582.4(0.5) 


7.9(0.4) 


0.58(0.08) 


6580.3(2.3) 


12.7(4.2) 


0.76(0.23) 


6582.5(0.9) 


7.6(1.3) 


0.27(0.07) 


6585.3(1.1) 


5.8(1.1) 


0.16(0.06) 










BA 








6553.5(5.4) 


215(30.0) 


8.22(3.2) 


6526.4(15.5) 


132.0(16.7) 


3.59(0.5) 


6486.5(1.6) 


90.1(13.1) 


1.77(0.38) 
















NA 








6572.1(1.8) 


6.4(2.0) 


0.19(0.07) 


6574.9(2.3) 


9.7(3.5) 


0.61(0.35) 


6575.0(0.9) 


12.3(1.5) 


0.51(0.08) 


6579.8(5.8) 


9.1(3.5) 


0.14(0.06) 








He 1 6678 


NE 


6696.7(1.3) 


8.7(3.0) 


0.08(0.02 
































He 1 7065 + CI 7100 


BE 


7101.3(5.2) 


129.9(18.2) 


2.78(0.52) 
































He 1 7065 


NE 








7086.0(0.7) 


10.0(1.6) 


0.14(0.02) 


























He 1 7065 


BA 


























6935.0(4.6) 


110.1(10.6) 


1.07(0.10) 








Call 729177324 


BE 


























7326.3(2.6) 


67.8(14.3) 


0.41(0.13) 


7313.8(4.5) 


110.6(10.2) 


0.73(0.11) 




BA 


























7250.7(2.9) 


26.4(6.1) 


0.09(0.03) 








017773 + S17725 


BE 


























7806.8(10.4) 


87.3(9.9) 


0.84(0.22) 


7796.8(12.5) 


86.6(15.5) 


0.17(0.04) 


He 11 8236(2) 


BE 


8240.6(7.8) 


104.2(18.8) 


17.20(7.10) 
































OI8446+CaII8498* 


BE 


























8494.9(3.3)/ 




14.6(1.1) 


8522.1(25.6)/ 




8.74(0.90) 


Call 8542, Call 8662 




























8674.5(3.3) 






8664.7(30.6) 






019261 +SI9222 


BE 


























9252.4(2.1) 


84.7(10.6) 


1.32(0.30) 









NOT!:.- Values lor the IAVUM a;id 11u.\ arc dependant on both the continuum placeman and ik dcbkidmii proteinic' cspcaalK m eases of kau line blending. The errors are 
given in brackets. 

" Line shape, "BE"=Broad Emission, "KE"=Narrow Emission. "BA"=Broad Absorption. "NA"=Narrow Absorption 
h in units of 10" !i ergs cm" 2 sec"'. The absorption line fluxes indicate the dux "subtracted" I'rom the continuum 
c This is a blend of He 11 4686 and C III 4648/ N HI 4640 
,l This is a blend of He 1 5876, Ba II 5874 and Na I D lines 
f This is a biend of He 1 5876, Ba II 5874 and Na 1 D lines 

? This profile is complex and can not be fitted with a single Gaussian. So we measured the flux by integrating the spectra directly. Two narrow "peaks" are superimposed on Ihe 
broad emission profile at 6581.0(3.5) Aand 6606.7(3.5) A 

s This profile has two clear separate peaks at the wavelengths quoted above. The line flux is measured for Ihe whole feature. Separation of the different components of the Ca II 
IR triplet (Ca II 8498, 8542 , 8662) and of 1 8446 was not performed due to the close blend of Ihese lines 
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Table 5. Narrow lines observed at high resolution. All observations were on days 17.4 or 36.3 unless indicated otherwise. 
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r vv nivi 


Flux c 






(kni / s) 


(km /s) 




firm /c^ 


(km /s) 




Emission Lines 














[Nell I] 3868.75 








JQ/L^fl 1 "\ 






H8+He I 3889.843| 








_l_Qfifi/'l \ 






He I 3964.729 










72(20) 


98ffi 1 9"i 


He(3970.072) 








TOU ( ^llj 
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\QKQ(7\ 




fl &A(C\ 1A} 


H 7 (4340.47) 




H\o) 


u. i y ^u. 10 ) 








[O III] 4363.209 










Ft^f 1 9 s ! 


n 9fi/'n ^7^ 


He I 4471.502 


ool lo ) 


OU^OO ) 


V.ZZyU. 1U ) 






U.DOl^U. 10) 


[Fe III] 4658.05 














H/3(4861.330) 


ooo\o ) 


A r K(Q\ 


n Qfif'n 91 ^ 






1 ^AYn 91 \ 




840(5) 


OV{t ) 








99 fn 0^^ 


[O III] 5006.843 


OOD^O ) 


*±o\v J 


U.OOI^U. LL ) 


1 Q/i 
-f-041 ( J 


oy^y ) 




He I 5015.678 




QQ/1 q\ 

OTjylO J 


n i n/fi ^^^^ 
u. iy i^u.uo ) 


-hoOo^ 1 J 


O 1 l ID j 


1 n^/'n 9^\ 


[Fe III] 5270.40 








TOUO ^ i-'-t J 


60(15) 




[N II] 5754.59 


- 


- 


- 


+848(5) 


45(3) 


0.30(0.04) 


He I 5875.661 


854(5) 


40(3) 


0.52(0.07) 


+863(5) 


42(5) 


0.79(0.17) 


Ha (6562.799) 


OOOy.) ) 


oo\o J 


9 1 ^{f\ 




OD{( ) 




Ha e (6562.799) 








TOUO^U J 


ovjyo j 


1 s&fn 1 7^ 


[N II] 6583.45 










QC\( Aft} 


f\Q(0 f\A~) 


He I 6678.152 


8^1 




n AAfn 1 ^ 








He I 7065.25 


000\0 J 


50(5) 


n d°Jc\ run 








He I 10830-f 










no 


1 ^(\(D 1^ 
±.o\)yvj.oo j 


He I 108309 










<100 


u.y^u.ou y 


Absorption 


lines 














H8+He I 3889.843 








AA (7\ 


9^91 


1 1 1 f'n n9^ 


He I 3964.729 








0A(A\ 




n 31 f'n fi9\ 


He(3970.072) 








ZOyO ) 






He I 4026.209 








98^8^ 
— Zo\o ) 


97/' £ ?\ 




H<5(4101.74) 










ir\( 1 

OU^-LO j 


u.yu^u.uo j 


H 7 (4340.47) 








-22(5) 


36(3) 


1.13(0.07) 


He I 4471.502 


-28(8) 


13(10) 


0.79(0.10) 


-23(10) 


25(10) 


0.27(0.05) 


H/3(4861.33) 


-28(5) 


25(6) 


0.35(0.11) 


-26(5) 


42(4) 


1.78(0.25) 


He I 5015.678 


-25(4) 


16(3) 


0.21(0.05) 


-21(4) 


16(3) 


0.87(0.07) 


He I 5875.661 


-27(4) 


15(3) 


0.26(0.04) 


-25(5) 


19(3) 


0.61(0.04) 


H«J6562.799) 


-26(5) 


32(5) 


1.19(0.10) 


-33(5) 


45(14) 


1.06(0.44) 


Hafl 










-27(8) 


40(10) 


0.87(0.07) 


He I 10830 


1 








-100(60) 


<100 


2.28(0.55) 


He I 10830 


! 








-100(60) 


<100 


2.1(0.55) 



a The rest wavelengths were obtained from the NIST Atomic Spectra Database and the CRC Handbook of Chemistry and Physics, 79th 
Edition. 

b For emission lines, total redshift is given. For absorption lines, shift is with respect to SN centre of mass (+846.9 km/s). Errors are 
given in brackets 

c in units of 10 — 14 ergs cm - 2 sec — . The absorption line fluxes indicate the flux "subtracted" from the continuum. Errors are given in 
brackets 

d We adopted A=3888.652 A for this blend 
e Observed at 38.5 d 
I Observed at 44.6 d 
9 Observed at 55.6 d 



discussed below are with respect to our adopted SN centre 
of mass velocity of +847 km/s. 

The following general conclusions can be reached. 
Where the S/N was sufficiently high, we can say that the 
peaks and troughs of the H and He lines exhibited similar 



velocity shifts with respect to the SN centre of mass frame. 
On day 17.4, their emission peaks were at +8±3 km/s, mov- 
ing to +15+3.5 km/s on day 36.3. The red wings of the H 
lines appear to shift to the red by around 15 km/s between 
the 2 epochs. Similar behaviour may be apparent in some 
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Table 4. Line identifications for the infrared spectra of SN 1998S. 



Table 4 - continued 



DAY ID 


A„„„ t (A) 


A„„„z. (A) 


LS a 


FWHM 


FLUX 6 






(Emis.) 


(Absorp.) 




(A) 




2 8 




19470(20) 




BE 


220(55) 


0.85(0.25) 




Br a 


19494(6) 




NE 


29(8) 


0.1(0.05) 




Br 7 


21775(35) 




BE 


380(90) 


0.3(0.1) 




Br 7 


21735(6) 




NE 


45(15) 


0.06(0.03) 


4.8 


He I 10830+P 7 


10910(15) 




BE 


241(35) 


13.9(3.1) 




p 


10966(8) 




NE 


32 (8) 


0.8(0.3) 




He I 10830 


10869(8) 




NE 


32 (8) 


1.3(0.2) 




He II 11626? 


11653(10) 




BE 


110(25) 


4.0(0.8) 




P „ 


12846(20) 




BE 


153(42) 


6.1(1.5) 






12847(8) 




NE 


38 (8) 


1 4(0 4) 


6.8 


He T 1 0830+P 


10930(20) 




BE 


210(22) 


14.9(1.9) 




lit _L lUOJU 


10875(5) 




NE 


29 (4) 


1.8(0.3) 


8.9 




19469(20) 




BE 


133(40) 


1.1(0.3) 




Br 7 


21726(20) 




BE 


115(25) 


0.5(0.1) 


11.9 


He T 1 0830 

1_J_C J. 1UOJU 




10850(5) 


NA 


26(7) 


2.1(0.6) 




He I 10830 


10874(5) 




NE 


50(10) 


3.1(0.9) 




He I 10830+P 7 


10950(20) 




BE 


126(20) 


6.9(0.6) 




p 


10973(4) 




NE 


17(5) 


0.4(0.2) 




p , 


12868(10) 




BE 


101(20) 


3.2(0.5) 




P j 


12864(4) 




NE 


26(7) 


0.7(0.3) 


12.8 


Brx 


19470(30) 




BE 


96(25) 


0.5(0.2) 




Br 5 


19497(12) 




NE 


23(8) 


0.10(0.04) 




Br 7 


21770(25) 




BE 


225(65) 


0.35(0.15) 




Br 7 


21729(6) 




NE 


35(15) 


0.13(0.04) 


17.7 


He T 1 0830 




10856(4) 


NA 


27(9) 


2.1(0.8) 




p 


10975(4) 




NE 


23(9) 


0.7(0.4) 




jne i luoou 


10872(2) 




NE 


25(6) 


1.7(0.4) 




He T 10830+P 


10945(20) 




BE 


170(40) 


7.5(2.5) 




P r, 


12855(5) 




NE 


21(8) 


0.5(0.2) 




P , 


12858(25) 




BE 


150(50) 


3.0(1.2) 


44.6 


He T 1 0830+P 




10642(30) BA 


226(35) 


11.0(2.5) 




Hp T 1 0830+P 


10915(30) 




BE 


210(35) 


30.5(5.5) 




Hp T 1 0830 




10852(5) 


NA 


23(0.5) 


0.25(0.08) 




O T 1 1 287 /1 1 294 

V_y jl 1 IlO 1 / J_ _L c/rr 


11300(15) 




BE 


150(25) 


5.8(0.7) 




+S I 11306 














Si II 12047 


11935(11) 




BE 


280(15) 


14.5(1.9) 




XU„ T 1 1 SOS 














p. 


12817(10) 




BE 


216(15) 


17.9(1.5) 




Mo- T 1 5031 


15027(15) 




BE 


105(30) 


2.1(0.4) 






18196(26) 




BE 


300(45) 


2.7(0.4) 




p 

r a 


18763(20) 




BE 


341(30) 


16.7(1.5) 




p 

1 a 


18749(4) 




NE 


34(6) 


0.65(0.15) 




cr^ 


19460(20) 




BE 


285(40) 


3.5(0.6) 




Br 7 


21690l'25) 




BE 


335(40) 


3.3(0.4) 


55.5 


P,s 


1 00451'! 0"! 




BE 


205(35) 


14 9(5 0) 




P«5 




10044(4) 


NA 


25(9) 


1.7(0.6) 




He I 10830+P 7 




10600(20) BA 


130(30) 


8.5(1.8) 




He I 10830 




10851(5) 


NA 


17(8) 


1.4(0.5) 




He I 10830+P 7 


10910(20) 




BE 


205(35) 


41.2 (1.9) 




O I 11287/11294 


11300(10) 




BE 


165(35) 


10.0(2.5) 




+S I 11306 














Si II 12047 


11950(20) 




BE 


345(40) 


18.0(1.7) 




+Mg I 11828 














P/3 


12825(15) 




BE 


205(25) 


29.5(2.5) 



DAY ID 



A peafc (A) X peak (A) LS^FWHM FLU? 
(Emis.) (Absorp.) (A) 



108.6- 
126.6 



He I 10830+P 7 
O I 11287 
+S I 11306 
Si II 12047 
+Mg I 11828 

Mg I 15031 

Br series, 

Si I 15888, [Fe II] 

[Si I], Mg I 

P Q 

Ca I 19309-19961 

Br 7 

CO 



10910(15) 
11327(15) 

12017(37) 

12875(15) 
15053(20) 
BLEND 



18833(40) 
19585(65) 
21770(45) 
~23400 



BE 275(25) 35.0(7.0) 
BE 250(30) 14.5(3.5) 

BE 330(55) 6.5(2.0) 



BE 
BE 
BE 



300(30) 20(2.5) 
208(35) 2.3(0.4) 
9. 8(2.5)£] 



BE 430(50) 33.5(6.5) 

BE 480(90) 6.7(2.5) 

BE 490(85) 3.9(1.2) 

BE - >30 



NOTE. -Values for the FWHM and flux are dependant on both the 
continuum placement and the deblending procedure especially in 
cases of heavy line blending. The errors are given in brackets. 



a Line shape, "BE"=Broad Emission, "NE"=Narrow Emission, 
"BA"=Broad Absorption, "NA"=Narrow Absorption 

4 ergs cm - 2 sec — 1 . The absorption line fluxes 
subtracted" from the continuum 
This is the integrated flux from 15500-17700 A 



b in units of 10" 
indicate the flux 








-12000 -9000 -6000 -3000 3000 6000 0000 12001 

Velocity (km/s) 

Figure 9. Evolution of the Ha line profile. Also shown are the 
Gaussian components that were used to fit the observed line pro- 
files. The km/s velocity corresponds to the adopted +847km/s 
velocity for the SN center of mass (see Section 3.2) 
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1.2 
1.05 






-12000 -9000 -6000 -3000 3000 6000 9000 12001 

Velocity (km/s) 

Figure 10. Evolution of the He I 10830 line profile. Also shown 
are the Gaussian components that were used to fit the observed 
line profiles. The km/s velocity corresponds to the adopted 
+847km/s velocity for the SN center of mass (see Section 3.2) 



of the He I lines also. In contrast, the absorption troughs 
(which also showed no difference in velocity shift between 
H and He) remained at a mean shift of -25±3 km/s be- 
tween the two epochs. The blue extremities of the H lines 
extended to ~-55 km/s on both epochs, whereas the blue 
limit of the He I lines was ~-40 km/s. The red extremities 
of the H and He I lines lay at ~+60 km/s apart from Ha 
which reached about +90 km/s. Figure [ill also reveals that 
the absorption component in most of the H and He I lines 
became significantly deeper, while the emission component 
showed a much smaller intensity variation. Only in Ha did 
the depth of the absorption remain about the same. (This 
may also apply to He I 4471 A but the S/N is very low.) On 
day 17.4, the only detectable forbidden lines were those of 
[O III] and only [O III] 5007 A was of good S/N. The blue 
edge of this line lay at ~-40 km/s, similar to that of the 
He I lines. However the red edge extended no further than 
+ 15 km/s. By day 36.3, the blue edge remained at -40 km/s, 
but the red edge had moved to ~+40 km/s. As mentioned 
above, we suggest that this is due to light travel time effects. 
By this epoch, lines of [Nil], [Ne III] and [Fe III] were also 
measurable and these all showed about the same absolute 
redshift and width as the [O III] lines. The [O III] profiles 
showed little variation in peak flux between the two epochs. 

In summary, the red wings of most of the emission lines 
showed a movement to the red of around +10 — \-2Q km/s 
between the two epochs, while the blue limits of the forbid- 
den line emission profiles and the allowed line absorption 
profiles remained roughly stationary i.e. the extent of the 
lines effectively increased by +10 — 1-20 km/s between the 



1 1 1 1 1 1 1 




III III 


Ha(6562.799) 




III 1 


II III 


III II 


I He 1 587+661 

1 1 v III 




Velocity (km/s) 

Figure 11. Narrow line profiles observed in the 17. 4d (fainter 
line) and 36.3d (darker line) high-resolution spectra. The vertical 
dashed line corresponds to the adopted redshift for the SN centre 
of mass, 847+3 km/s (cf. Section 3.2). To help the comparison of 
the line profiles, the 36.3d continua have been displaced vertically 
to coincide with the 17. 4d continua. 



two epochs due the redshift of the red wing. We attribute 
this to light travel time effects. We also note that the fact 
that the [O III] 5007 A profile varies in time, together with 
the similarity of its shape and redshift on day 36.6 to those 
of other forbidden lines, indicates that the SN CSM is in- 
deed the origin of these lines. This rules out the alternative 
line-of-sight emission line region suggested by Bowen et al. 
(2000). 

The narrow P Cygni lines of He I 10830 A acquired 
at UKIRT (resolution ~100 km/s) on days 44 and 55 are 
also listed in Table |^. The emission profiles and peak-trough 
profiles are unresolved. The emission peak lies at velocities 
of +29+20 km/s on day 44 and -7+15 km/s on day 55, 
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which is the same as for the optical He I lines, to within 
the uncertainties. The line strengths are virtually identical 
at the two epochs. While the main emission and absorption 
components are unresolved, there is some indication that 
the blue wing of the absorption is barely resolved, reaching 
a blueshift of ~-300 km/s i.e. similar to that found for the 
wing of the broad Ha absorption. 

We conclude that the narrow lines arose in a CSM wind 
of velocity of 40-60 km/s. Highest velocities (~60 km/s) are 
seen in the H lines, and lowest (~40 km/s) in the forbidden 
lines. The greater velocity extent of the H lines may be in 
part due to the stronger effects of thermal broadening in 
this element. For a temperature of 10,000 K, the thermal 
broadening contribution has a FWHM of 21 km/s in hydro- 
gen compared with just 5 km/s in oxygen. It is also possible 
that radiative acceleration effects could contribute to the 
greater width of the allowed lines. Acceleration of the CSM 
could be driven by the UV radiative precursor and/or cos- 
mic ray precursor (Fransson, Lundqvist & Chevalier 1996), 
with the innermost part of the CSM being accelerated most. 
The optical depth in allowed lines in the accelerated CSM 
component could be large (due to high radiative excitation) 
thereby increasing the width of the scattered lines. However, 
it is quite plausible that the volume emission measure of the 
accelerated component is low and so would not contribute 
to the forbidden lines. 

Turning briefly to the broader absorption component in 
the optical H and He I lines, we can see (Figure [ll]) that this 
feature strengthened somewhat between days 17.4 and 36.3. 
Similar behaviour is reported for the H lines by Bowen et 
al. (2000). However, Bowen et al. also note the presence of a 
similarly broad absorption in Mg II 2796, 2803 A which (a) 
may have declined in strength between days 32 and 39, and 
(b) never exhibited narrow P Cygni lines. We also note that 
the Ha and H/3 lines are sufficiently well defined to allow 
us to conclude that the broad minima were blueshifted by 
~-115 km/s on day 17.4, moving redward to ~-85 km/s by 
day 36.3. 



4 DISCUSSION 

4.1 The circumstellar medium of SN 1998S 

The CSM of SN 1998S has been described and discussed by 
a number of authors (LOO, GOO, Fassia et al. 2000, Bowen 
et al. 2000, Lentz et al. 2000). In the following discussion 
we adopt the general consensus about the form of the CSM 
and the evolution of the spectra. In particular, we confirm 
the two-phase CSM proposed by LOO, but we deduce some- 
what different, and occasionally firmer constraints on the 
size and density of the CSM. In this paper we focus mostly 
on what can be deduced from the first-season observations. 
Later observations will be described in a subsequent paper 
(Fassia et al., in preparation). We envisage that before it ex- 
ploded, the progenitor star of SN 1998S had undergone two 
phases of mass loss, with a gap or a phase of reduced mass 
loss in between. This gave rise to two shells of CSM. Follow- 
ing LOO, we shall refer to these as the Inner CSM (ICSM) 
and Outer CSM (OCSM). The interaction of the ejecta with 
these two shells is responsible for the early (t <10days) 
and late (t >40days) time appearance of broad emission 



lines lacking a P Cygni absorption component (cf. Chugai 
& Danziger 1994, Chevalier & Fransson 1994). 

4-1.1 The inner circumstellar material (ICSM). 

We recall that the early-time broad component had a 
FWHM>4000 km/s, was present as early as day 3, and 
had faded significantly by day 12-14 (Qiu et al. (1998); this 
work). While lines of this width in SNe are usually asso- 
ciated with the ejecta, the absence of any corresponding 
P Cygni absorption is a robust indicator that these lines 
are driven by the high energy radiation of the shock re- 
sulting from the interaction of the ejecta with dense CSM 
in the immediate vicinity of the supernova (Chugai 1990). 
The ratio of the luminosities of the Ha and H/3 broad com- 
ponents in the first week spectra is /_ff a //if^=1.45±0.25 
and the ratio of the fluxes of the P/3 and H7 lines is 
/p/3//h7=0.22±0.07. If we take into account the redden- 
ing assuming Av = 0.68j^ 35 (Fassia et al. 2000) then these 
ratio become /h c «//h^=L1±0.7 1 /p l a//H7=0.1±0.06. These 
low ratios indicate a high degree of thermalisation for these 
lines and thus imply a high density and a large optical depth 
for the line-forming region. 

The fading of the broad emission components by 
days 12-14.3 and the appearance of broad absorption com- 
ponents in H and He I by days 14.3-17.7 (cf. Figure is 
strong evidence that by that time the ejecta had completely 
overrun the ICSM. The highest ejecta velocity was about 
~10,000 km/s relative to the SN center of mass. We there- 
fore estimate that the outer limit of the ICSM was <90 AU. 
By ~15 days the photospheric radius was already ~75 AU 
(Fassia et al. 2000). 

4.1.2 The region between the ICSM and OCSM. 

As can be seen in Figures |?], ^|the spectra during 14-27 days 
are dominated by broad P Cygni- type absorption profiles. 
During this phase, the extreme blue edge of the broad ab- 
sorption indicated a velocity of -8,000 km/s relative to the 
SN centre of mass. We adopt this as the maximum velocity 
of the ejecta. The ejecta collided with the OCSM sometime 
after day 27. The high resolution spectra of day 36 show no 
sign of broad emission. However, by day 44, strong broad 
emission had appeared in the IR lines. We deduce that the 
8,000 km/s shock reached the inner edge of the OCSM on 
day 40±4 indicating a radius of 185±20 AU. So the region 
between the ICSM and OCSM extended from <90 AU to 
~185 AU. 

There is evidence that even during the phase when the 
ejecta had overrun the ICSM but not yet reached the OCSM, 
the ejecta/CSM interaction contributed significantly to the 
spectrum. While the SN 1998S spectral lines at this phase 
are similar to those observed in the spectra of normal SNe II 
during their photospheric phase (eg. Blanton et al. 1995, 
Fassia it et al. 1998), the SN 1998S spectrum differs in that 
it has a much stronger, bluer continuum and shallower ab- 
sorption profiles. As mentioned in Section 3.1.3, Branch et 
al. (2000), Lentz et al. (2000) and LOO suggest that the low 
contrast of the broad lines relative to the continuum is due 
to external illumination of the line-forming region by light 
from the CSM-ejecta interaction. We agree with this sug- 
gestion and furthermore propose that the observed strong 
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continuum emission arises in a massive, cool, dense shell of 
ejecta formed at the ICSM-ejecta discontinuity (e.g. Cheva- 
lier & Fransson 1994). Such a shell is likely to be produced 
by the interaction of the ejecta with the ICSM, and would 
have a velocity jump of ~ 10 3 km s^ 1 between the inner 
part of the ejecta and the shell. This could have led to en- 
hanced continuum emission in two ways. In the first mech- 
anism, even after the ejecta have completely overrun the 
ICSM, the existence of the velocity jump would have re- 
sulted in a reverse shock wave which ran through the ejecta 
and which would have remained powerful enough to gener- 
ate strong continuum emission from the dense shell. For the 
second mechanism we suppose that between the ICSM and 
the OCSM there actually existed a tenuous, uniformly dis- 
tributed circumstellar medium which we shall to refer to as 
the MCSM (Mid-CSM). Interaction of the dense, swept- up 
shell with the MCSM would have released power (kinetic 
luminosity) according to L = 2-kR 2 pv A or L oc R 2 oc t 2 . In- 
terestingly, between days 17.7 and 44.6 the P/3 flux grew ~6 
times, in good agreement with the expansion factor between 
these two epochs viz. (t'2/ti) 2 — 6.3. We note, however, that 
if the P/3 lines arose from the dense shell and were thermal- 
ized (saturated) at a constant excitation temperature then 
the expansion of the dense shell alone could account for the 
L oc R 2 behaviour. Nevertheless, if the dense shell did ex- 
ist then it may have contributed to the broad line emis- 
sion which appeared in the IR by day 44 i.e. at this epoch 
the line luminosities may have been composed of contribu- 
tions from both radioactive and ejecta/CSM energy sources, 
where "CSM" includes ICSM, MCSM and OCSM. However, 
determination of the relative luminosities from these sources, 
and what fraction of the ejecta/CSM luminosity arose from 
the different zones, is beyond the scope of this paper. 

4.. 1.3 The outer circumstellar material (OCSM). 

The presence of the OCSM is revealed in two phases. The 
earlier is due to its interaction with the UV/X-ray flash from 
the SN, and the later is due to its interaction with the SN 
ejecta. 

a) The early phase. 
We attribute the early, narrow, usually unresolved, lines to 
recombination and heating following ionisation of the slowly- 
moving OCSM by the UV/X-ray flash at shock break-out. 
These lines were present in the earliest spectra (day 3.3) as 
the narrow (FWHM<500 km/s) components of the emission 
lines. By days 12-14, the narrow H, He I lines began to take 
on the form of P Cygni profiles, and this is clearly confirmed 
in the day 17.4 WHT/UES spectra. However, the high res- 
olution spectra also show that the narrow lines in H and 
He I are actually made up of at least two components. The 
components are best described as a P Cygni-like profile of 
maximum velocity ~50-60 km/s and a broader, blueshifted 
absorption extending to around -350 km/s. The high reso- 
lution spectra also show narrow forbidden lines extending to 
±40 km/s of the centre of mass velocity. 

We consider first the narrowest lines or line components 
i.e. those extending to ±40-60 km/s. We envisage that the 
narrow P Cygni-like profiles in H and He I were produced by 
two effects. The recombination cascade in the OCSM would 
generate narrow emission lines. In addition, the recombina- 



tion would populate excited levels in the OCSM, resulting 
in photospheric photons being resonantly scattered, hence 
forming classical P Cygni lines. Thus, the narrow line emis- 
sion component comprises photons emitted both through 
recombination and resonant scattering. We note that the 
He I 6678, 7065 A lines appear only in emission and only 
in the 17.4 day spectra. He I 6678 A is produced by a sin- 
glet transition to the 2p 1 Po level. The lack of an absorption 
in this line is presumably due to rapid decay from 2p 1 Po 
to the metastable 2s 1 S level. In contrast, the 5015 A line 
is produced by a transition directly to the 2s 1 S level, and 
that accounts for the absorption in this line. Thus, He I 
6678 A is a virtually pure recombination line and its dis- 
appearance by day 36.3 is presumably due to a lowering of 
the recombination rate by then. He I 7065 A is produced 
by a triplet transition to the 2p 3 Po level. The presence of 
absorptions in the 5876 and 4471 A lines (which both also 
decay to the 2p 3 Po level) indicates a significant population 
in 2p Po. However, the cross-section for resonant scatter- 
ing from this level via the 7065 A transition is only about 
10% and 50% relative to that of the 5876 A and 4471 A 
transitions respectively. Given that the 4471 A line barely 
exhibits an absorption component above the noise, it is not 
surprising that there is little evidence of absorption in the 
7065 A line. The stronger pronounced emission component 
in 7065 A relative to 4471 A may be due to its being of lower 
excitation energy. Collisional excitation from the metastable 
2s 3 S level would thus be expected to produce a significantly 
stronger 7065 A line. The disappearance of the 7065 A line 
by day 36.3 is probably due to cooling of the electron gas. The 
forbidden lines would have been excited as a result of the 
heating of the OCSM electron gas following recombination. 
The appearance of narrow [N II] , [Ne III] and [Fe III] lines 
between days 17.4 and 36.3 can be attributed to the for- 
mation of the relevant species following recombination from 
higher ionisation. 

We can estimate the CSM density in the region where 
the [O III] lines were formed by using the ratio of their fluxes 
R=(/a4959 ± /a5007)//a4363- L00 measured the ratio of the 
narrow [O III] lines present in their 4d and 5d spectra days 
and found that for temperatures in the range 15000 K<T 
< 30000 K the electron density of the CSM was 6 x 10 6 < 
n c < 2 x 10 7 cm" 3 . However, the SPSDAS package of IRAF 
that they use to derive these values is appropriate only for a 
low density (n< 10 5 cm -3 ) plasma. In addition, theoretical 
models developed to study the temperature and structure of 
the CSM around supernovae (Lundqvist & Fransson 1988) 
have shown that the CSM temperature is very sensitive to 
the ionising spectrum of the supernova and varies both with 
radius and in time. The temperature range considered by 
LOO for the O Ill-rich gas is much narrower than in the 
models of Lundqvist & Fransson (1988), which indicate 10 4 - 
10 5 K for a distance greater than ~200 AU in a situation 
similar to that of SN 1998S. The upper limit, 10 5 K, is set by 
the fact that O III is collisionally ionised to O IV at higher 
temperatures. 

To calculate the CSM density from the [O III] line ratio 
we therefore assumed that 10 4 K < T < 10 5 K and used the 
6- level [O III] model atom described in Maran et al. (2000). 
This model includes O III atomic data from Aggarwal (1993) 
and Galavis et al. (1997) and takes into account collisional 
de-excitation that plays an important role for n > 10 5 cm~ 3 . 
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Figure 12. Solutions of the 6-level model atom for values of the 
[O III] line ratio, R, estimated from the SN 1998S spectra (cf 
Section 4.1.3) 

Unfortunately, our 3.3 d-17.7 d low-resolution spectra are 
of insufficient S/N for us to estimate the ratio, R, of the 
oxygen lines. However, we were able to estimate the ratio 
on day 36.3 using our high resolution observations. Before 
de-reddening, the ratio is R=3.3 ± 1.0. De-reddening using 
the Cardelli et al. (1989) law with A v = 0.68±£-jg (Fassia et 
al. 2000), changes the ratio to R=2.9±1.5. Therefore, using 
the 6-level model atom mentioned above we estimate that 
for 1.4<R<4.4 the lower limit to the CSM density is 1 x 10 6 
cm" 3 . A CSM of this density and velocity (40-50 km/s) 
is consistent with it being a red supergiant wind (Dupree 
1986, Salasnich, Bressan & Chiosi 1999, Jura & Kleinmann 
1990). No solution for the density is found for T < 8000 K 
(see Fig. 12), which means that this is a firm lower limit to 
the temperature in the O Ill-rich gas at this epoch. If the 
temperature is T < 11000 K, which cannot be excluded from 
the results of Lundqvist & Fransson (1988), no upper limit 
to the density can be given. The day 36.3 line ratio therefore 
only places lower limits on the density and temperature of 
the O Ill-rich gas. 

For 17.4 d, only an upper limit is available for the flux 
of [O III] 4363 A from our high resolution spectra and so 
we get R> 1.9 (dereddened). From Figure |l^ we see that 
this, unfortunately, does not constrain the density for the 
temperature range 10 4 -10 5 K. 

We turn now to the broader (~350 km/s) shallower ab- 
sorption component seen in the H and He I lines. We re- 
call that this feature deepened between days 17.4 and 36.3. 
The minima also moved redward from ~— 1 15 km/s to ~- 
85 km/s. The discrete difference in its FWHM, relative to 
that of the narrow component might indicate that the two 
components were formed in different regions. The presence 
of the broader absorption as late as day 35 implies that it too 
must have arisen in the OCSM. It could not have been due 
to the ICSM since it had been totally overrun by then. In ad- 
dition, an origin in the region between the ICSM and OCSM 



seems unlikely since not only must the density of material 
there have been quite low (see above) but also by day 35 the 
undisturbed region between the shock and the inner edge of 
the OCSM was very thin - barely ~25 AU. While a wind ve- 
locity ~350 km/s is too high for a red giant, blue supergiant 
(BSC) winds do exhibit velocities in the range ~100 km/s 
to ~3000 km/s (eg. Abbott et al 1981). Thus a possible ex- 
planation for the broader absorption is that it occurred in a 
component of the OCSM produced when the progenitor was 
going through a BSC phase. The deepening and shift to the 
red suggests the effect of a BSC wind "piling up" against 
the inner boundary of the RSG wind. This would produce 
a negative velocity gradient in the BSC wind, as well as 
making its density gradient less steep. As the SN shock pro- 
gressed, it would first engulf the faster, more tenuous part of 
the BSC wind, leaving the outer, slower-moving zone undis- 
turbed, thus accounting for the decrease in the absorption 
blueshift. The deepening absorption could be due to an in- 
crease in the density of the outer, slower-moving zone of the 
BSC wind as it became increasingly compressed between 
the shock and the RSG wind boundary. The division of the 
OCSM into two zones in this scenario would also provide a 
natural explanation for the apparently discrete difference in 
the FWHM of the narrower and broader absorption features. 

A possible alternative explanation for the ~350 km/s 
absorption component is that it might be produced by the 
acceleration of the innermost OCSM by the radiation pres- 
sure of the UV precursor, mentioned at the end of Sec- 
tion 3.3. This would also produce the required negative ve- 
locity gradient. Details of the physics of the high velocity 
wind component will be considered in a separate paper. 

b) The later phase 
As indicated above, we take the presence of strong, broad 
Paschen and Brackett lines on day 44 to imply that the 
shock/OCSM interaction had already already begun by this 
date. This is further supported by the weakening of the nar- 
row components of the H and He I lines (see Figs. 9 and 
10), which we attribute to the invasion of the OCSM by 
the shock. The broad lines persisted to at least day 692 
(Fassia et al. in preparation) with material still moving as 
fast as 4600 km/s relative to the SN centre of mass at that 
time. From this we infer that the OCSM extended at least 
as far as 1800 AU from the supernova. For a wind veloc- 
ity of 50 km/s, this implies that the the ejection of the 
OCSM began at least 170 years ago, in agreement with 
LOO. We agree with LOO that line ratios at this later time 
continue to imply a high density in the line formation re- 
gion of the OCSM. Our estimate for the deredened [O III] 
ratio is (/ A 4959 + /a5007)//a4363 = 2.5 ± 0.9 (average of 
values from days 72 and 97 ) which suggests (see Fig- 
ure |l2[ ) densities n e > 1.5 x 10 6 cm -3 for a temperature 
of lCT K < T < 10 5 K, consistent with the density derived 
at earlier epochs by us and by LOO. 

From the above we can obtain lower limits for the 
OCSM mass and for the mass loss rate that produced it. 
Assuming that the OCSM has the form of a spherically sym- 
metric wind of constant velocity v w . the density would vary 
as oc r~ 2 . Consequently the mass of the OCSM is: 

M > 47vrlm H nfi(— - 1) (1) 
where fi is the mean atomic weight in amu, thh is the mass 
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of the hydrogen atom, r out , ri n are the inner and outer limits 
of the OCSM and n is the total number density at the in- 
ner limit. Making the conservative assumption of complete 
ionisation and that the OSCM consists mostly of hydro- 
gen, we have n ps 2n e and p=0.5. Therefore, substituting 
n e ~ 1.5 x 10 6 cm" 3 , r in ~ 185 AU and r ou t > 1800 AU, 
we have M > O.OO3M0. For a wind velocity of 50 km/s 
and a constant mass loss we obtain a mass loss rate of 
> 2 x 10 M©/yr taking place for at least 170 years be- 
fore the explosion. 

4.2 Carbon Monoxide 

Molecular emission from a supernova was first observed in 
the form of CO emission from SN 1987A (Catchpole & Glass 
1987; McGregor & Hyland 1987). Measurement of the first 
overtone spectrum as early as day 110 was reported by 
Meikle et al. (1989). However, it is possible that CO formed 
even earlier. The sudden rise in the K — M colour at about 
day 100 (Bouchet & Danziger 1993) may have been due to 
CO fundamental emission in the M-band. The evolution of 
the SN 1987A first overtone emission spectrum was moni- 
tored by Meikle et al. (1989, 1993) to 1.5 years. Other species 
responsible for molecular emission from SN 1987A include 
SiO together with more tentative identifications of CS, FeO, 

and HeH + (Dalgarno, Standi & Lepp, 1997, and refer- 
ences therein). Of these molecules, emission from CO has 
been most extensively studied. 

Analysis of the first overtone CO band in SN 1987A 
has demonstrated the value of this type of study in deter- 
mining the velocity, excitation temperature and CO mass in 
the CO-emitting zone. The first models (Spyromilio et al. 
1988), which assumed LTE, provided good fits for days 192 
& 255 implying a CO mass of ~10" 4 M Velocities of 
2000 km/s and 1800 km/s, respectively, were inferred for 
the two epochs. At later epochs poorer fits were obtained. 
This was attributed to deviation from a Boltzmann distribu- 
tion. Subsequent interpretations by several groups involved 
non-LTE modelling for the vibrational transitions. An ex- 
ample is the work of Liu & Dalgarno (1995) who obtained 
much-improved fits for the late-time data. They inferred a 
mass of CO as high as 0.04 Mq on day 110 but which was 
reduced quickly to -~3x 10 -3 M Q by day 250 due to destruc- 
tion by the fast electrons. They found that, when present, 
the vibrational transitions of CO provide the dominant cool- 
ing mechanism of the C/O zone during the first 2 years, 
reducing the temperature to ~700 K by the end of that 
period. Prior to SN 1998S, the only other supernova (be- 
sides SN 1987A) for which CO emission has been detected 
is the type II SN 1995ad at around 105 days, observed by 
Spyromilio & Leibundgut 1996. Following Liu & Dalgarno, 
they suggest that a CO mass as high as 0.05 Mq might 
have formed. In SN 1998S intense first overtone CO emis- 
sion spectra were measured at two epochs (days 110 and 240, 
with our definition of t = 1998 March 2.68 UT) by GOO and 
at a single epoch (day 109) by us. We shall compare their 
conclusions with ours at a later point in this section. 

4-2.1 Non-LTE model of CO emission. 

To interpret the CO spectrum of SN 1998S, we have con- 
structed a model following a similar non-LTE approach to 



that of (Liu & Dalgarno 1995), including CO formation and 
gamma-ray energy deposition. The only additional feature 
is a simple treatment of non-local radiative coupling of the 
first overtone lines in a clumpy envelope. By adjusting the 
model to reproduce our day 109 CO spectrum (Fig. g, |l3| ) 
we can estimate the mass, density, temperature and velocity 
of the CO and C/O-rich matter as well as the mass of 56 Ni. 
To compare the model with the data, the model luminosity 
was converted to flux assuming a distance of 17 Mpc (Tully, 
1998). 

In the model we assume that the CO emission arises 
from an ensemble of C/O-rich clumps embedded in a homo- 
geneous (on average), macroscopically- mixed core. All the 
C/O-rich clumps are composed of equal proportions by mass 
of C and O, and have similar size and physical conditions. 
Following Liu & Dalgarno (1995) we assume that the total 
mass of the clumped C/O-rich mixture is M c \ = O.2M0. 
This is reasonably consistent with mass estimates from stel- 
lar evolution computations for stars in the range 12-25 Mq 
(Heger 1999). Gamma-ray energy deposition is determined 
using the escape probability approximation with an absorp- 
tion coefficient of fc 7 = 0.03 cm 2 g _1 . However, it is likely 
that the 56 Ni is also clumped. Consequently, as the gamma- 
rays travel from the 56 Ni clumps to the C/O-rich clumps 
they will suffer attenuation, the amount of which will de- 
crease with age. This process has been invoked in SN 1987A 
to account for the slow evolution of the CO luminosity at 
early epochs (Liu & Dalgarno 1995). To allow for this ef- 
fect in SN 1998S, we follow Liu & Dalgarno by including 
an energy reduction factor tj) by which the deposition rate 
in the C/O-rich material is multiplied. To estimate ip for 
SN 1998S on day 109, we compared the first overtone CO 
flux with that reported by GOO on day 225. The density of 
the C/O-rich clumps (p c i) may be higher than the average 
density in the core (p c ), and so we specify the overdensity 
with a density contrast parameter \ = Pci/p c > 1. The CO 
concentration is primarily controlled by CO formation via 
the radiative association C + O — > CO + hv and destruction 
via ionization by fast electrons CO + e — > CO + + e followed 
by the dissociative recombination CO + + e — > C + O (Liu & 
Dalgarno 1995). 

The fraction of deposited energy channelled into 
Coulomb heating was obtained using calculations of Kozma 
& Fransson (1992). The electron temperature (T c ) of the 
C/O-rich clumps was found assuming energy balance be- 
tween the energy deposited as heat by Compton electrons 
and the radiation emitted in the CO fundamental and first 
overtone bands. Cooling in [OI] 6300, 6364 A lines was found 
to contribute typically less than 4% of that due to the CO 
bands, in agreement with the analysis by Liu & Dalgarno 
(1995). 

The non-LTE population of vibrational levels was calcu- 
lated assuming that radiative and collisional transitions be- 
tween neighbouring vibrational levels are dominant. This is 
a sensible approximation since the radiative transition rates 
for the first overtone are an order of magnitude smaller than 
those for the fundamental. The electron collision rate coef- 
ficients were taken from Liu et al. (1992), while oscillator 
strengths are from Kirby-Docken & Liu (1978). Twelve vi- 
brational states, with 101 rotational levels for each vibra- 
tional state, were considered. This yields a line list of 2000 
lines for the first overtone band and 2200 lines for the fun- 
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damental band. The CO radiative rates are affected by con- 
tinuum radiation originating outside the clumps. The effects 
of this were determined by assuming that the core luminos- 
ity is equal to the overall deposited radioactive luminosity. 
The external continuum spectrum was approximated by a 
dilute blackbody with a colour temperature of 5000 K. In 
practice, the model is not very sensitive to the choice of the 
temperature since the effect of the external continuum ra- 
diation is small compared to that of the collisions. This is 
in accordance with the Liu & Dalgarno (1995) conclusions 
about SN 1987A. After the populations of the vibrational 
levels were found, the rotational levels of individual vibra- 
tional states were populated according to the Boltzmann 
distribution. 

Radiative transfer in the clumpy CO matter was treated 
in terms of the average escape probability. This includes the 
combined effects of the escape probability of photons from 
an average CO cloud, together with absorption by the other 
CO clouds in the core. Our models reveal that variation of 
the CO cloud radius between 10 km s _1 and 100 km s _1 
(assuming homologous expansion) only slightly affects the 
shape of the spectrum. We adopted a cloud radius of a = 
50 km s _1 , very similar to the oxygen clump size estimates 
for SN 1987A (Stathakis et al. 1991). The weak sensitivity of 
the spectrum shape to the cloud size is due to the opposing 
trends of photon absorption in the CO emission source cloud 
and in other clouds. 

The final model is specified using four free or "input" 
parameters: core velocity (v c ), core mass (M c ), density con- 
trast (x), and 56 Ni mass (Mm)- It should be stressed that 
not all these parameters are equally constrained by the ob- 
servations. In particular, good fits over a wide core mass 
range may be achieved by simply adjusting the other three 
parameters appropriately. However, we find that models 
with a core mass M c < 0.4 Mq require an implausibly high 
density contrast in the C/O-rich matter (x > 10). On the 
other hand, models with a core mass M c > 4 Mq require a 
56 Ni mass > 0.2 Mq, inconsistent with the bolometric lumi- 
nosity of SN 1998S around day 130 (Fassia et al. 2000). We 
have therefore considered the two extreme cases set by these 
considerations viz. model A with a core mass of 0.4 Mq and 
model B with a core mass of 4 Mq . 

4-2.2 Comparison of the CO model with the day 109 
spectrum. 

We find that both models are able to reproduce the observed 
CO first overtone emission for plausible combinations of the 
free parameters v c> X an d Mm- From these models we ob- 
tain physically reasonable values for the "output" parame- 
ters Mco, T e , rib {rih is the density of the C/O-rich matter 
expressed as a baryon number density) . In Fig |l3| we show 
model A fits for core velocities of 1600, 2200 and 2800 km/s. 
We can immediately eliminate the two extreme values. In 
the 1600 km/s model, the visibility of the quasi-periodic 
structure produced by the R-branches of the 0-2 and 1-2 
vibrational transitions is too high. In the 2800 km/s model 
the visibility is too low and the rising blue edge too shal- 
low. We deduce that the velocity of the CO emission zone is 
2200 ±300 km s . A similar result is obtained for model B. 
For the optimum velocity of 2200 km/s, in Table ^ we show 
the input and output parameters for the two models. 
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Figure 13. Carbon monoxide first overtone model A (thick line) 
(cf. Table ^) compared with observations (thin line). Models are 
shown for three CO emission zone velocities: 1600 km s — 1 (upper 
panel), 2200 km s _1 (middle) and 2800 km s _1 (bottom). 



We find that rib, Mco and T e are quite similar for mod- 
els A and B. To narrow down further the range of models en- 
compassed by A and B, additional information is required. 
One such fact is the estimate by Fassia et al. (2000) of a 
56 Ni mass of 0.15 ± 0.05 Mq, based on the bolometric light 
curve of SN 1998S. This is consistent with model B but not 
model A (Table |^) . A caveat is the possible contribution to 
the luminosity from the ejecta-wind interaction which might 
decrease the 56 Ni mass implied by the bolometric light curve. 
However, Fassia et al. point out that the observed bolometric 
decline rate matches well the radioactive luminosity decay 
rate, arguing against a significant ejecta-wind contribution 
to the luminosity. Another argument in favour of model B 
is based on the issue of the energy reduction factor. The 
model A fit requires ip — 0.2 — 0.3, indicating a substan- 
tial difference between the average intensity of gamma-rays 
m the 56 Ni zone and C/O-rich clumps. Yet, in this model 
the optical depth of the core for gamma-rays is only about 
1. Moreover, macroscopic mixing should lead to even lower 
optical depths in the 'walls' between 56 Ni and C/O zones. 
These low optical depths are inconsistent with a ip value as 
low as 0.2-0.3, as is required to fit the data. On the other 
hand, the version of model B rm which fits the observed 
spectrum has a core optical depth of ~10 on day 109, con- 
sistent with its ip value of ~ 0.5 (Table ^J). 

In summary, we favour the high mass core model 
(model B) with a core mass of 4 Mq, core velocity of 
2200±300 km/s, electron temperature of 3650±200 K, neg- 
ligible overdensity of the C/O-rich clumps (density contrast 
parameter x ~1) an d a 56 Ni mass of about 0.2 Mq. (Note 
that the C/O is still clumped in the sense of chemical abun- 
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Table 6. Carbon monoxide model parameters 



Input parameters (see text) Output parameters (see text) 

Model M c M m X i> n b M C o T c 

(M ) (M Q ) (10 11 cm- 3 ) (10- 3 M Q ) (K) 

A 0.4 0.07 8 0.24 0.77 1.2 3700 

B 4 0.2 1 0.52 0.97 1.3 3650 



dance.) Unfortunately, the mass of the mixed core does not 
permit us to draw definite conclusions about the progenitor 
mass. The reason is that we lack information about the rel- 
ative proportions of metal-rich and H/He-rich materials in 
the core. That the amount of H/He-rich material might be 
significant is supported by the observation of sharp-topped 
Ha and He I 10830 A line profiles at t ~ 100 days, indicating 
the presence of low velocity hydrogen and helium (v ~ 1000 
km s" 1 ). 

GOO have analysed the CO emission from SN 1998S cov- 
ering about the same epoch as did our study. Matching their 
day 110 CO first overtone spectrum with a range of models 
they deduce a CO velocity of 2000-3000 km/s. They argue 
that a CO velocity of 3000 km/s would imply a progeni- 
tor mass in excess of 25 Mq, which is outside the plausible 
progenitor mass range for our adopted C/O-rich core mass. 
However, we believe that our analysis eliminates the upper 
half of their velocity range, which presumably would reduce 
the inferred progenitor mass. But in any event, as pointed 
out above, there are difficulties in reaching firm conclusions 
about the progenitor mass. We note that GOO's CO velocity 
value depends quite strongly on the presence in their spec- 
trum of a "smooth rise to the blue edge of the CO feature" . 
However, our spectrum, taken within a day of their's does 
not show this smooth rise. Instead, our spectrum remains 
roughly flat to about 2.27 /im and then rises abruptly. In 
Figure [m] we compare the two spectra, binned to 25 A per 
pixel. It can be seen that the agreement is excellent up to 

2.27 /jm, but then the spectra become more discrepant. In 
particular, the GOO spectrum exhibits a "bump" at about 

2.28 /mi which contributes significantly to the perception of 
their smooth rise to the blue edge. However, their spectrum 
appears to be rather noisy redward of 2.3 /im and it may 
be that noise is also responsible for their 2.28 /im bump. 
Our value for the electron temperature, 3450-3850 K, also 
confines the 3000-4000 K range suggested by GOO. 

As indicated above, Liu & Dalgarno (1995) deduce a 
CO mass of ~ 10" 2 M in SN 1987A at 110 days, which 
is an order of magnitude greater than for SN 1998S at a 
similar epoch (cf. Table 1). The reason for this is that the 
56 Ni mass and the energy reduction factor, ip, are both x3 
lower for SN 1987A than for SN 1998S, leading to ~ 10 more 
efficient CO formation in SN 1987A. 



5 CONCLUSIONS 

We have presented and discussed optical and infrared spec- 
tra of the type Iln SN 1998S, covering epochs from a few 
days after explosion to over 100 days later. Our observa- 
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Figure 14. Comparison of the first overtone CO spectrum ob- 
tained on day 109 (solid line) with the spectrum obtained by GOO 
on day 110 (dotted line). The spectra have not been corrected for 
the redshift of the supernova. 



tional coverage of SN 1998S makes a significant contribution 
to the study of the type Iln phenomenon. Of particular note 
is our acquisition of (a) contemporary spectra in both the 
optical and IR bands at a range of epochs, and (b) the most 
extensive set of high resolution spectra ever for this type of 
supernova event. The spectroscopic evolution of SN 1998S 
was complex. It can be understood in terms of the inter- 
action of the supernova with a two-component progenitor 
wind which we have referred to as the ICSM and OCSM. 

Collision of the ejecta with the ICSM accounts for the 
early spectral features. From the fading of the broad emis- 
sion components by ~14d we deduce that the outer bound- 
ary of this wind lay at less than 90 AU from the centre. Es- 
timation of the inner wind velocity is difficult due to (a) the 
ongoing interaction of the shock with this wind even at the 
earliest of times, and (b) the lack of high resolution spectra 
at these times. We can only say that the early-time spectra 
exhibited features of width less than 400 km/s. However, 
these may have arisen in the OCSM. If, like LOO, we assume 
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that the inner and outer winds had similar velocities then 
from our high resolution measurements of the OCSM we can 
infer that the ICSM wind commenced less than 9 years ago. 
However, the possibility of an RSG-BSG evolution would 
argue against this assumption. 

Examination of the spectra indicates that the OCSM 
extended from 185 AU to over 1800 AU. Our high resolu- 
tion spectra have revealed that the OCSM has a velocity of 
40-50 km/s. Assuming a constant velocity during this time, 
we can infer that it commenced more than 170 years ago, 
and ceased about 20 years ago. During this period the out- 
flow was high - at least 2 x 1O _5 M0 yr _1 , corresponding to 
a mass loss of at least 0.003 Mq. An outflow of this strength 
and velocity is similar to those seen in cool supergiants. The 
broader absorption feature (^350 km/s) in H and He I may 
have arisen from a component of the outer CSM shell pro- 
duced when the progenitor was going through a later blue 
supergiant phase. Alternatively, it may have been due to 
UV-precursor radiative acceleration of the inner part of the 
OCSM. 

Our analysis of the CO emission together with the bolo- 
metric light curve Fassia et al. 2000) also indicates a mas- 
sive progenitor, with a mixed metal/He core of M ~ 4 Mq 
i.e. comparable to that of SN 1987A. SN 1998S is only the 
third core-collapse SN for which post- 100 day AT-band spec- 
tra have been published, and yet all SN three events have 
exhibited first-overtone CO emission. (Fassia et al. (1998) 
also published IR spectra for the core-collapse SN 1995V. 
No CO was identified but the latest K-h&nd spectrum was 
only at 69 days.) A picture is therefore beginning to emerge 
where CO plays a ubiquitous role in the evolution of core- 
collapse SNe at late times. In particular, the powerful cooling 
property of CO is likely to lead to conditions in which dust 
may condense in the ejecta. Excess IR emission has been re- 
ported in SN 1998S on days 225, 260 and 255 (GOO), day 253 
(Garnavich et al. 1998) and for a series of dates to as late 
as day 691 (Fassia et al. 2000), and this may be due to dust 
condensation. (An IR excess at 130 days was also reported 
by Fassia et al. (2000) but they argue that it was unlikely 
to have arisen from dust condensation). In a future paper 
(Fassia et al. in preparation) we shall examine in detail the 
evolution of the late-time IR excess. 

LOO used spectropolarimetry at early times to infer sig- 
nificant asymmetry in the OCSM and in the ejecta-ICSM 
interface. While our discussion of the high resolution OCSM 
lines assumed spherical symmetry, we have not provided 
quantitative explanations for the relative shifts between the 
allowed and forbidden lines. It may be that this will also 
require the introduction of an asymmetric model. GOO also 
deduced asymmetry in the OCSM, on the basis of the broad 
line profiles of H and He on days 240, 275, 370. We shall 
discuss these profiles in a future paper. 
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